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Abstract. Methods of optimal control are applied to a model system of interacting 
two-level particles (e.g., spin-half atomic nuclei or electrons or two- level atoms) to 
produce high-fidelity quantum gates while simultaneously negating the detrimental 
effect of decoherence. One set of particles functions as the quantum information 
processor, whose evolution is controlled by a time-dependent external field. The 
other particles are not directly controlled and serve as an effective environment, 
coupling to which is the source of decoherence. The control objective is to generate 
target one- and two-qubit unitary gates in the presence of strong environmentally- 
induced decoherence and under physically motivated restrictions on the control field. 
The quantum-gate fidelity, expressed in terms of a novel state-independent distance 
measure, is maximized with respect to the control field using combined genetic and 
gradient algorithms. The resulting high-fidelity gates demonstrate the feasibility of 
precisely guiding the quantum evolution via optimal control, even when the system 
complexity is exacerbated by environmental coupling. It is found that the gate duration 
has an important effect on the control mechanism and resulting fidelity. An analysis of 
the sensitivity of the gate performance to random variations in the system parameters 
reveals a significant degree of robustness attained by the optimal control solutions. 
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1. Introduction 

The transfer of information between elements of a quantum computational system 
requires the use of entangling quantum interactions Undesired interactions between 
the system and its surroundings can destroy quantum coherences and thus are a 
critical obstacle to successful quantum computation (QC). The feasibility of creating 
high-fidelity quantum gates in the presence of environmentally-induced decoherence is 
one of the most important problems to overcome for practical QC. In particular, in 
spin-based solid-state realizations of QC [2J [3j HI [5] one encounters a difficult task 
of effectively separating a multiparticle quantum system into interacting and non- 
interacting components. 

Quantum error correction (QEC) enables fault-tolerant QC [6], but only when 
the errors in quantum gate operations are sufficiently small [7]. Therefore, it is very 
important to decrease the errors caused by decoherence. This problem has inspired 
significant interest in various methods of decoherence management, including the use 
of decoherence-free subspaces and noiseless subsystems [H [9j QUI [HI H2L quantum 
dynamical decoupling [T31 [141 H3 IlBl HZ], [TBI [19l [20], schemes based on stochastic control 
[2Tj . optimal control techniques [221 E31 EH ES, EE! [27] , and multilevel encoding of logical 
states [28] . 

The method of optimal control [291 130] enables managing the dynamics of complex 
quantum systems in a very precise and specific manner and therefore is especially 
useful in QC. In addition to applications to the problem of dynamical suppression of 
decoherence [221 E31 El E3 ESI E7], optimal control theory (OCT) [3U [32] was also 
successfully used to design unitary quantum gates in closed systems [331 EH EH [361 [37]. 
The optimal control of quantum gates in the presence of decoherence still remains to be 
fully explored. In [28] we previously considered the optimal control of quantum gates for 
qubits encoded in multilevel subspaces; this method makes quantum gates immune to 
mixing and decoherence that occur within the encoding subspaces. Recent works [381139] 
developed specific techniques, involving optimizations over sets of controls operating in 
pre-designed "weak-decoherence" subspaces. In the present paper we propose a different 
approach in which the full power of OCT is used to generate the target gate with the 
highest possible fidelity while simultaneously suppressing strong decoherence induced 
by coupling to a multiparticle environment. This method does not rely on any special 
pre-design of the system parameters to avoid or weaken decoherence (e.g., using multiple 
levels as in [28j, tunable inter-qubit couplings as in [38], or auxiliary qubits as in [39]); 
the only control used in the present approach is a time-dependent external field. 

A similar OCT-based approach was recently used [Hj] to design quantum gates for 
solid-state qubits in the presence of decoherence. However, the objective in [ID] was to 
optimize a purity- dependent quality factor (or, in [38], the purity itself), instead of the 
actual gate fidelity. In the present work we demonstrate that although improving the 
purity of the quantum information processor (QIP) is necessary for performing a high- 
fidelity quantum gate, it is not sufficient. Even if the QIP is completely decoupled from 
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its environment at a given time, this does not ensure that the desired gate operation 
will be performed at the decoupling time. Therefore, we optimize a gate fidelity [JT] 
which directly measures the distance between the target quantum gate of the QIP and 
the actual transformation of the composite system. Optimization techniques were also 
applied recently to QEC [321 S3]- In contrast to QEC, our approach does not require 
ancilla qubits and is not limited to the weak decoherence regime. The optimal control 
of quantum gates can potentially be used in conjunction with QEC to achieve fault 
tolerance with an improved threshold. 

In this work, we consider a model system composed of interacting two-level particles, 
for example, spin-half atomic nuclei or electrons or two-level atoms. A small set of 
particles serve as qubits in the QIP; the rest of the particles serve as an effective 
environment. The qubits are directly controlled by a time-dependent external field, while 
the environmental particles do not directly couple to the field. The control objective is to 
generate target quantum gates in the QIP with the highest possible fidelity. The optimal 
control field must perform the desired gate operation while simultaneously suppressing 
the qubit-environment interaction and restoring lost coherence to the QIP. This model is 
sufficiently simple to allow for a full numerical treatment of the entire composite system, 
and the results are relevant to important physical applications, in particular, to spin- 
based solid-state realizations of quantum gates [2j El SI E] . For example, our model bears 
a similarity to systems in which an electron spin (or a pair of electron spins) is coupled 
to a nuclear spin bath [51 HH H5] . Coherent manipulation of electron spins via rapid 
electrical control of the exchange interaction has been successfully demonstrated in such 
systems [5]. The analysis reported in the present work indicates that the employment 
of the optimal control methods may increase the effectiveness of coherent management 
of coupled spin dynamics. 

The paper is organized as follows. Section [2] presents the model (including an 
explicit matrix form for the simplest case of one qubit coupled to a one-particle 
environment) and schemes of multiparticle couplings. In section[3j we consider a distance 
measure that quantifies the fidelity of quantum gates. This fidelity is independent of 
the initial state and is evaluated directly from the evolution operator of the composite 
system. Section H] investigates the dynamics of decoherence in the uncontrolled system 
for various values of system parameters. In order to fully explore the utility of OCT, we 
select a set of parameters that enhances the loss of coherence in the uncontrolled system. 
In section EJ we describe in detail the genetic and gradient optimization algorithms. 
The results obtained with the optimal controls are presented and discussed in section |6j 
Section [7] investigates the robustness of optimal solutions to uncertainties in the system 
parameters. Finally, section [8] concludes with a summary of the results and discusses 
future directions. 
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2. The model system 

We use a model of iV interacting two-level particles (e.g., spin-half particles or two- 
level atoms), which are divided into the QIP, composed of m qubits, and an n-particle 
environment (N = m + n). The qubits are directly coupled to a time-dependent 
external control field, while the environment is not directly controlled and is managed 
only through its interaction with the qubits. The evolution of the composite system 
of qubits and environment is treated in an exact quantum-mechanical manner, without 
either approximating the dynamics by a master equation or using a perturbative analysis 
based on the weak coupling assumption. The Hamiltonian for the composite controlled 
system, H = H + H c + H int , has the form (h = 1) 

N m N-l N 

h = uiS iz - iiiC(t)s ix - 7^ • s^. (i) 

i=l i=l i=l j>i 

Here, Sj = (Si X , Si y , Si Z ) is the spin operator for the ith particle (Sj = \<y%-, in terms of the 
Pauli matrices), H is the sum over the free Hamiltonians LUiSi z for all iV particles (a;, is 
the transition angular frequency for the ith particle), He specifies the coupling between 
the m qubits and the time-dependent control field C(t) are the dipole moments), 
and H mt represents the Heisenberg exchange interaction between the particles (7^ is the 
coupling constant for the ith and j'th particles). This model is particularly relevant to 
spin-based solid-state realizations of quantum gates [21 [31 HI [5] . 

Now consider the simplest case of one qubit and a one-particle environment 
(m = n = 1) in more detail. The Hamiltonian in this case is: 

H = uiS lz + c^SW - iiC{t)S lx - 7S1 ■ S 2 , (2) 

where 7 = 712. We use the orthonormal basis: 

|1> = |+> 2 , |2) = |+) 1 ®|-) 2> |3) = |-)i®|+) 2 , |4) = |-)!®|-) 2 , (3) 

where Si Z \±)- t = ±||±)j. The Hamiltonian (T5]) in the basis ([3]) has the following matrix 
form: 



2 



( wi + uj 2 - §7 -nC(t) \ 

LOi - uj 2 + §7 -7 -/iC(t) 

-fj,C(t) -7 u 2 - LOl + 57 

^ -fiC{t) -ui - uj 2 - §7 J 



(4) 



In addition to the simplest case of a two-particle system described above, we also 
consider situations where one qubit is coupled to a multiparticle environment (m = 1 
and n — 2, 4, 6). For m — 1, the coupling constants are given by 

= f 7, for i = 1 and j = 2, . . . , JV, ( , 

llj 1 0, for 2 < i < N, U 
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which means that the qubit interacts with each environmental particle with the same 
coupling constant 7, and the environmental particles are not directly coupled to each 
other. For n = 2, the system can be modeled by a linear chain with the qubit q± at the 
center, equally coupled to both environmental particles e 2 and e 3 : 

e 2 < — > qi < — > e 3 (6) 

For n = 4, the system can be modeled by a two-dimensional lattice with the qubit q± at 
the center, equally coupled to four environmental particles {e 2 , . . . , e 5 }: 

e 4 
I 

e 2 < — ► qi < — ► e 3 (7) 
I 

Similarly, for n = 6, the system can be modeled by a three-dimensional lattice with 
the qubit at the center, coupled to six environmental particles. In these lattices, it is 
assumed that the Heisenberg interactions decay exponentially with distance [2], and 
therefore environmental particles on the vertices of the square (n = 4) and cube (n = 6) 
are neglected. 

A different model with nearest-neighbor couplings is also considered in the case of 
n — 4. The system is modeled by a linear chain of particles, with the qubit at the center 
and each particle coupled only to its nearest neighbors with the same coupling constant 
7: 

e 4 < — > e 2 < — > qi < — > e 3 < — > e 5 (8) 

The case where two qubits are coupled to a one-particle environment (m = 2 and 
n — 1) is used to develop an entangling quantum gate (specifically, the controlled-NOT 
gate) in the presence of a simple environment. This system can be modeled by the 
following two-dimensional triangular lattice: 

e 3 

713/ \723 ( 9) 

712 

qi < — ► q2 

where the two qubits are denoted as q\ and q 2 , and the environmental particle as e 3 . 
Such a model is relevant, for example, for a dilute nuclear spin bath [3]. Values for this 
set of coupling constants are given in section 16.41 

3. The distance measure 

Our objective is to generate an evolution of the QIP which at some time tf will be as 
close as possible to the target quantum gate. The problem of evaluating the actual gate 
fidelity is complicated by the fact that the evolution of the QIP is non-unitary due to the 



Optimal control of quantum gates and suppression of decoherence 



6 



interaction with the environment. Nevertheless, it is possible to define a useful measure 
of the distance between the target quantum gate of the QIP and the actual evolution 
operator of the composite system [JT]. 

Let U(t) G XJ(2 N ) be the unitary time-evolution operator of the composite system 
and G G U(2 m ) be the unitary target transformation for the quantum gate of the QIP 
(where U(cf) denotes the group of all d x d unitary matrices). The evolution of the 
composite system is governed by the Schrodinger equation, 

U(t) = -\H{t)U{t), (10) 

with the initial condition U(0) = I 2 n (where Id denotes the d x d identity matrix). The 
gate fidelity depends on the distance between the actual evolution U = U(tf) at the final 
time tf and the target transformation G. In order to perform a perfect gate, it suffices 
for the time-evolution operator at t = tf to be in a tensor-product form U opt = G <E> $, 
where $ G U(2 ra ) is an arbitrary unitary transformation acting on the environment ||| 
Therefore, the following objective functional is proposed [JT] as the measure of the 
distance between U and G: 

J = \ N mm{\\U-G® $|| | $ G U(2 n )}, (11) 

where || • || is a matrix norm on the space (C) of d x d complex matrices (in the 
present case d = 2 N ), A at is a normalization factor, and J is minimized over the set 
of all unitary $. It is useful to expand G, and U in orthonormal bases. Let {\i)}, 
{|^)} and {\i) (g) \u)} be orthonormal bases that span the Hilbert spaces of the QIP, 
environment, and composite system, respectively. The corresponding expansions read 

2 m 2 n 

G = $ = E $ -'IW(^l> (12a) 

i,i'=l u,u'=l 
2 m 2" 



U= Y] V U u ,\i){i'\ ® \v)(i/\. (12b) 
• •/ i / 1 vv ' 

1,1 =1 v,v — 1 

Using in ( fTTl) the Frobenius norm, defined as 

||X|| Fr = [Tr(XtX)] 1/2 VXGM d (C), (13) 



and A at = 2 ( iV + 1 )/ 2 ; the distance measure becomes 

J = 

where Q G M 2 ™(C) is given by 

2" / 2" 



1/2 



1-2-^Tr (Vgtg) , (14) 



f,i/'=l \i,i'=l 



| We do not consider in the present work a more general situation where the composite system itself 
is open and may not be unitary. 
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Since < J < 1, it is convenient to define the gate fidelity as F = 1 — J. An important 
property of this distance measure is its independence of the initial state. In contrast to 
some other distance measuresj§| J is evaluated directly from the evolution operator U, 
with no need to specify the initial state of the system. This property of J reflects our 
objective of generating a specified target transformation for whatever initial state, pure 
or mixed, direct-product or entangled. 

Note that in the ideal case when there is no coupling to the environment, i.e., the 
QIP is a closed system with unitary dynamics, the distance measure (tHJ) becomes 

J= [l-2- m |Tr(Gt[/ q )|] 1/2 ) (16) 

where U q = U q (t{) is the unitary evolution operator of the QIP at the final time. Another 
distance measure used in the literature [37] for closed systems is J cs = 1— 2~ m \Tr(G^U q ) |, 
i.e., J cs = J 2 . For example, in section [6] we report optimization results which, in the 
case of closed QIP systems, are J ~ 10~ 6 and J ~ 10~ 4 for one- and two-qubit gates, 
respectively, corresponding to the values J cs ~ 1CT 12 and J cs ~ 10~ 8 , respectively. 

4. Decoherence dynamics of the uncontrolled system 

The loss of coherence in the QIP, caused by the interaction with the environment, 
is detrimental to the quantum gate performance. In order to better understand the 
mechanism of optimal control, we first study the decoherence process in the uncontrolled 
system. The state of the QIP at time t is described by the reduced density matrix: 

p q (t) = Tr cnv [p(t)} , (17) 

where p{t) is the density matrix of the composite system and Tr env denotes the trace 
over the environment. A useful measure of decoherence is the von Neumann entropy 

Eh 

S vN (t) = -Tr{p q (t)ln[p q (t)}}. (18) 

For a pure state, SVn = 0, while for a maximally mixed state of a /c-level system, 
SVn = ln(/c). We explore the decoherence dynamics of the QIP by studying the time 
evolution of the entropy S v ^(t) for the uncontrolled system (in this section) and under 
the influence of optimal time-dependent control fields (in subsequent sections). The 
initial state used for the entropy calculations is 

m N 

l*o) = (g)|->i® (g) \+h ( 19 ) 

i=l j=m+l 

(i.e., initially all qubits are in the state |— ) and all environmental particles are in the 
state |+)). Recall that the distance measure J of (TT4"|) is independent of the initial state 

§ Relationships between various distance measures, including some presented in |46j and generalizations 
of (fT4|) . are discussed in more detail in [41] . 
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and consequently so are the optimal control fields found for the target gates and the 
corresponding fidelities. We choose some initial state only for the entropy calculations, 
which are done to illustrate the decoherence dynamics after the time-evolution operator 
is determined (for either a controlled or uncontrolled system). Therefore, the specific 
choice of the initial state ( f!9l) places no limitations whatsoever on the generality of the 
optimal control results. 

We set the unit of time, thereby introducing a natural system of units, by arbitrarily 
choosing u>i = 1 for all simulations (this implies that one period of the first qubit's free 
evolution is 2tt). Details of the dynamics depend on the system parameters (i.e., the 
frequencies and coupling constants for the uncontrolled system). In the simplest case of 
the uncontrolled system of one qubit coupled to a one-particle environment (m — n— 1), 
the initial state is |\&o) — | — )i ® |+)2, and the time evolution can be solved analytically: 



= e" i7t/4 | cos (fit) | — h> +isin(Ot) 
p q (t) = cos 2 (fit) |-)(-| + sin 2 (fit) 



2fi 1 1 2fi' 1 



(20) 



(21) 



4fi 2 1 " 1 4fi 2 

where we use a simplified notation: | — h) = |— )i ® |+)2, |H — ) = |+)i <8> |— )2, 
and fi = |[(o>i — oj 2 ) 2 + 7 2 ] 1//2 is the Rabi frequency. Due to discreteness of the 
environment's spectrum, the loss of coherence is reversible. If the transition frequencies 
are degenerate, U\ = u 2 , then the state of the composite system, |\I/(t)), oscillates 
between two direct-product states, |— +) and \-\ — ). In this case, complete coherence 
revivals will occur whenever sin (fit) = or cos (fit) = 0, i.e., at times tf eg) = /br/(2fi) 
{k G N). However, if ui\ ^ u 2 , then \^{t)) oscillates between the initial direct- 
product state | — h) and an entangled state (a superposition of | — h) and | H — )). 
Therefore, complete coherence revivals will occur only when sin(fit) = 0, i.e., at times 
tk = kn/Q (k G N). If \ui — u 2 \ <C 7, then, in addition to the complete revivals 
at times t^, partial revivals will occur at times t ( k part) « (k - |)7r/fi (k G N). The 
maximum loss of coherence depends on the values of 7 and \ui — o^l- For a given value 
of 7, closer frequencies enhance the interaction between the qubit and environment, 
causing higher peak values of decoherence (i.e., the entropy) and longer revival times. 
Figure [1] shows the time-evolution of the entropy for the uncontrolled system of one 
qubit and a one-particle environment, with 7 = 0.02, iO\ = 1, and various values of 
uj 2 . The entropy dynamics shown in figured], obtained by numerically propagating the 
Schrodinger equation fflOl) . and are in full agreement with the analytical results above. 
In particular, we find the first-revival times t\ ~ {50.0, 140.7, 313.2} for uj 2 = (jr — x) _1 
and ti w {43.9, 136.1, 313.2} for u 2 = tt - x with x = {2, 2.1, 2.14}, respectively. These 
values fully agree with the analytical formula for tf. obtained above. Also, for x = 2.14, 
the frequency difference \ui — uj 2 \ ~ 0.00159 is about one order of magnitude smaller 
than 7, and, correspondingly, a partial revival is found numerically at t[ part ' ) w 156.6, in 
agreement with the analytical result. 

For the optimal control simulations below, the system parameters are chosen to 
ensure complex dynamics and strong decoherence: values ofj/u are up to 0.02, which 
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is significant for QC applications, and the frequencies uJi are close (but not equal), 
to enhance the interaction. For one qubit coupled to a one-particle environment 
(m = n = 1), we choose 

lo 1 = 1, uj 2 = (tt- 2. 14)- 1 « 0.99841. (22) 

Imposing upper limits on the gate duration (tf < 60) and coupling constant (7 < 0.02) 
places the dynamics of the uncontrolled system in the regime where decoherence 
increases monotonically with time (before the entropy reaches its maximum value of 
SVn ~ In 2). This dynamical regime approximates some of the effects that the QIP 
would experience from a larger environment, in particular, preventing restoration of 
coherence to the qubit by uncontrolled revivals. Thus, any increase in coherence may 
be attributed exclusively to the action of the control field. 

When selecting the parameters of a multiparticle environment, we apply the same 
criteria for maximizing decoherence of the uncontrolled system, as described above. 
Figure [2] illustrates the uncontrolled time-evolution of the entropy for a one qubit coupled 
to n-particle environments (n = 2,4,6), with 7 = 0.02. The frequencies of the qubit 
and pairs of the environmental particles are given by 

oj\ = 1, (23a) 
ujj = (tt — Xj)^ 1 , = n — Xj, j — 2,4, ... , n, (23b) 

{2.14, n = 2, 

2.14,2.1, n = 4, (23c) 
2.14,2.1,2, n = 6. 

For example, for n = 4, the frequencies of the four environmental particles are 
approximately {0.96007, 0.99841, 1.00159, 1.04159}. 



5. Optimal control algorithms 



In the context of optimal control, the objective is to maximize the fidelity of the target 
quantum gate over a set of time-dependent control fields. The target quantum gates 
considered in this paper include the Hadamard (H t ), identity (I2), phase (vr/8), and 
controlled-NOT (CNOT) transformations: 




CNOT 




1 

exp(i7r/4) 



(24a) 



(24b) 



Collectively, H t , ir/8, and CNOT constitute a universal set of quantum gates for QC 
PQ. Identity is included to preserve an arbitrary quantum state during a specified time 
interval, e.g., while operations are performed on other qubits. 
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In maximizing the gate fidelity, we employ a combination of two optimization 
techniques, a genetic algorithm and a gradient algorithm. For a given target gate, the 
genetic algorithm first locates a parameterized control field that achieves a reasonable 
value of fidelity (e.g., F > 0.95), then the gradient algorithm further improves this result 
by lifting the parameterization restriction on the field. This section describes the details 
of these search algorithms. 



5.1. Optimization with the genetic algorithm 

When the genetic algorithm is used, the gate fidelity F is maximized with respect to a 
parameterized control field 

m 

C(t) = f(t)^2 AiCos(uit + 9i), 0<t<t { . (25) 
i=i 

Here, f(t) is an envelope function incorporating the field's spectral width, tf is the 
gate duration, and A+, Ui, and 9i are the amplitude, central angular frequency, and 
relative phase of the zth component of the field, respectively. A combination of these 
optimization parameters (called "genes") represents an "individual" whose "fitness" is 
defined as the fidelity of the gate generated by the corresponding field. A collection of 
individuals constitutes a "population" (we use population sizes of ~ 250). At each 
generation, we evaluate the fitness of all population members and create the next 
generation by crossover and mutation of genes of the fittest individuals (crossover 
and mutation rates are between 20 and 40 percent). A novelty of this algorithm 
implementation is the inclusion of the control duration tf as one of the optimization 
parameters. 



5.2. Optimization with the gradient algorithm 

Removing the constraints on the control field imposed by the parameterized form (1251) 
provides the potential for more effective control of the system. In this case optimal 
control fields are found by minimizing the following functional 



K 



J + Re f ' Tr { \(l(t) + iH(t)U(tj\ B{t))dt+^ f ' |C(t)| 2 dt. 
Jo J ' 2 J 



(26) 



In addition to the distance measure J of (114j) . K includes a constraining term and 
a cost term. Upon minimization of K, the first integral constrains U(t) to obey the 
Schrodinger equation (B(t) is an operator Lagrange multiplier) and the second integral 
term penalizes the field fluence, 



£ 



\C{t)\ 2 dt, 



(27) 



with a weight a > 0. 
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5.2.1. Optimal control equations. An optimal control field is obtained by solving a 
set of equations that follow from the variational analysis of K as a functional of B(t) 
and U(t). Here, we derive the corresponding functional derivatives [IE] and boundary 
conditions. The functional derivative of K with respect to B(t) yields 



SK 



Re 



U(t) + \H(t)U(t) 



1 T 



(2f 



so that the condition 5K/5B(t) = results in the Schrodinger equation ( TTTJT) for U{t). 
Next we compute the functional derivative of K with respect to U(t): 



5K 
5U(t) 



Re 



f 5J 



+ B 1 (t t )5(t - t t ) - B{t) - iB{t)H{t) 



T 



(29) 



Since J depends only on U = U(t{), we obtain 5J/5U(t) = (dJ/dU)5(t — t{). Therefore, 
the condition 5K/ 5U (t) = results in two equations: 



B(t) = 
B T (t{ 



iB(t)H(t), 

_dJ_ 
' ~ ~dU' 

We will also use the functional derivative of K with respect to C(t), 

5K 



5C(t) 



Im {Tr \fiU(t)B(t)]} + aC(t), 

rn 



(30) 
(31) 

(32a) 
(32b) 



to guide the gradient search, as described in section 15.2.21 below. 

The initial condition for U(t) is U(0) = I 2 n and the final condition for B(t) is given 
by (13T1) . In order to find the explicit form of dJ/dll, first consider a scalar function 
y(Z(x)), where Z is a matrix function of the scalar variable x. Using the chain rule, we 
obtain 

dy dZl K ^JdydZ T ' 



dy_y- dy dZ KK 
dx / -— ' dZ KK i dx 



E 



dZ KK i dx 



Tr 



dZ dx 



Setting y = Tr(Z), implies that 



Now let Z = (Q f Q) V2 and x = U, 



dy = 
dZ 



(33) 



(34) 



a b (a complex scalar variable). The matrix indices a 

1/2 

and b range from 1 to 2 N . Note that (Q^Q) is not an analytic function of U a b, but 
it can be expressed as an analytic function of U a b and U* b . Therefore, a generalized 
complex derivative [48J is applied to calculate dZ/dx, so that U* b and subsequently 
are treated as constants when differentiating (Q^Q) 
find that 

dZ _ d(Qtg) 1 / 2 
dx 



dU, 



ab 



l/ 2 



-1/2 



with respect to U a b- Thus we 
dQ 



Q 



dU, 



(35) 



ab 
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By combining (I3"3"j) - (13"5]) . we obtain 



dy 



dx dU, 



Tr 



ab 



tnW 2 



{Q j Q) 



Tr 



WQ)- 1/2 Q'— 



dU, 



ab 



(36) 



With the above notation, J = (1 - 2- N y) 1 / 2 . Noting that (dJ/dU) ab = dJ/dU ab and 
using (136|) . we finally derive 



fd£ 



-N 



Tr 



{Q ] Q) 



'1/2 



Tr 



dQ 
dU, 



G 



ab 



fa/2"] \b/2 n ] 



ni\a mod 2 n )(b mod 2 n |. 



(37a) 
(37b) 



Equation (137bj) is obtained from (fT5l) . using the fact that A; mod k = k. In (137bl) . the 
states are elements of the environment's orthonormal basis {|z^)}, and \x] denotes the 
smallest integer greater than or equal to x. The explicit form of the boundary condition 
for B(t) is obtained by substituting f|3T|) into (|3ip . 



5.2.2. The numerical procedure. Optimal control fields are found using an iterative 
gradient algorithm described below. An initial guess for the control field is needed at 
the first iteration. Typically, we use the output of the genetic algorithm as the initial 
guess for faster convergence, although fields of the form (125]) with a random choice 
of parameters can be used as well. At each iteration, U(t) is propagating forward in 
time with the Schrodinger equation ( ITUi) and the initial condition U(0) = I 2 n. The 
resulting matrix U = U(tf) is used to determine the final condition (13T!) for B(tf). Then 
B{t) is propagated backward in time with the time-reversed Schrodinger equation ( |30l) . 
All propagations are performed using a toolkit for computational efficiency [19]. The 
resulting U(t) and B(t) are utilized to compute the functional derivative SK/SC(t) of 
( 132]) . which then adjusts the control field for the next iteration. The adjustment of the 
control field for the kth iteration (k G N) is given by 



C«(t) = C (fc - 1} (t) - psm r (wt/tf) 



5K 



5C(t) 



c(t)=c( k - 1 Ht) 



(38) 



where < (3 < 1 and | < r < 1 are constants used to modify the magnitude of the 
field adjustment. The multiplier sin r (nt/t{) ensures that the control field C(t) is nearly 
zero at the initial and final time, which is a reasonable physical restriction on the field. 
This iteration routine continues until we observe no further improvement in K, which 
manifests the achievement of convergence. 

Despite the lack of direct coupling of the control field to the environment, it can 
be shown that the composite system described by ([T|) is completely controllable (up to 
a global phase), as defined in [50J. However, the restrictions on the gate duration and 
on the shape of the control field limit the achievable fidelity. 
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6. Results of optimal control in the presence of decoherence 

6.1. One qubit coupled to a one-particle environment 

We consider the optimally controlled Hadamard, identity, and phase gates generated for 
a single qubit coupled to a one-particle environment (m = n = 1). Fidelities for these 
one-qubit gates are presented in figure [3] for various values of the coupling constant 7. 
The control fields optimized for the actual values of 7 result in fidelities above 0.9991. 
In particular, for the Hadamard transform, we obtain F > 1 — 10~ 6 for 7 = (a closed 
system) and F ~ 0.9995 for 7 = 0.02 (the strongest coupling considered). In contrast, 
when the control field optimized for 7 = is applied to the system with 7 = 0.02, it 
generates a gate with a poor fidelity, F ~ 0.9063. This result demonstrates that optimal 
solutions designed for the ideal case of a closed system have little value when applied 
to realistic open systems. However, the optimal control algorithm is able to generate 
quantum gates with very high fidelities, if coupling to the environment is explicitly taken 
into account. 

The optimal control fields that generate the one-qubit gates (with a one-particle 
environment and 7 = 0.02) are shown in figure HI These fields are intense, with 
maximum amplitudes larger than 2.0 (in the units of h — u)\ — fa — 1). The gate 
duration is tf = 25.0 (about four periods of free evolution). The exact time structure 
of an optimal field is not intuitive and is tailored to the particular control application. 
For example, control fields optimized for 7 = 0.02 are not only more intense than those 
optimized for 7 = 0, they also have very different structures. One common feature of 
the optimal control fields presented in figure H] is that they are approximately symmetric 
about t ~ tf/2. We suggest that this property of the fields is related to the reversibility 
of the system dynamics: the periods in which the information flows from the QIP to the 
environment are followed by periods in which the information flow is reversed, in order 
to restore the coherence of the QIP. 

Figure [5] shows the time behavior of the von Neumann entropy of the QIP for 
optimally controlled one-qubit gates (with tf = 25.0 and 7 = 0.02). By comparing 
figures [5] and [IJ we observe that the optimal control dramatically enhances coherence 
of the qubit system in comparison to the uncontrolled dynamics. Decoherence is 
suppressed by the control at all times, but especially at the end of the gate operation 
(i.e., for t = tf). For example, S^fff) < 10~ 7 for the Hadamard gate with 7 = 0.02, 
which means that at t = t{ the qubit system and environment are almost completely 
uncoupled. Inspecting eigenvalues of the controlled Hamiltonian, we find that the intense 
control field creates significant dynamic shifts of the energy levels. Specifically, under 
the influence of the optimal control field, four of the six transition frequencies of the 
composite system experience high-amplitude oscillations (following the corresponding 
changes in the field strength). This effect is mainly responsible for reducing the qubit- 
environment interaction during the control pulse. However, achieving extremely low 
final-time entropies and correspondingly high gate fidelities requires the employment of 
an induced coherence revival. For the selected set of the system parameters, revivals in 
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the uncontrolled dynamics occur at times much longer than if (specifically, i[ part ^ pa 156.6 
and if ~ 313.2), so that the almost complete coherence revival observed at i = if is 
induced exclusively by the control field. 

For very short gate durations (if < 5), a different type of optimal solution is found. 
The control fails to induce revivals at such short times and therefore generates gates 
with smaller fidelities (e.g., F ~ 0.9874 for the Hadamard transform with 7 = 0.02 and 
if m 2.33). In this short-time regime the control relies on the decoherence suppression 
via dynamic shifting of the energy levels and on very fast operation (trying to perform 
the target transformation in the shortest time possible to limit the effect of decoherence), 
but not on the creation of coherence revivals. Such short-time controls can be useful 
for environments with very dense spectra, for which the induced-revival times will be 
impractically long. 

We study in detail how the choice of the control duration if affects properties 
of the optimal control field, gate fidelity, and decoherence dynamics. Specifically, we 
optimize the one-qubit Hadamard gate (with a one-particle environment and 7 = 0.02) 
for all integer values of if between 2 and 40 using the gradient algorithm described in 
section 15.21 For if < 5 we find the fast-control no-revival regime described above. 
Interestingly, most optimal control fields with if > 5, in addition to inducing an 
almost complete coherence revival at the final time, also produce a partial revival at 
approximately if/2. Optimal control fields with 5 < if < 20 typically exhibit large 
amplitudes and fluences and strong low-frequency components. For if = 25 we find 
the optimal control field that generates the quantum gate with a better fidelity while 
having a smaller amplitude and fluence, as compared to the fields obtained for shorter 
control durations. As if increases to 25, the gate fidelity increases to approximately 
0.9995, the final-time entropy decreases to approximately 10~ 7 , and the maximum field 
amplitude decreases to approximately 2.0. However, increasing if above 25 does not 
improve the optimal gate performance; the field amplitudes, gate fidelities, and final- 
time entropy values change very slightly for 25 < if < 40. The physical interpretation 
of this behavior is that the control requires some time (if > 25 in the present case) to 
almost completely reverse the information flow between the QIP and environment, and 
induce a nearly perfect coherence revival. From these results, it appears that the pulse 
duration is a very important characteristic of the control fields employed for quantum 
gate generation. 

6.2. The Kraus-map dynamics of the qubit 

The time-dependent state of the QIP, which is coupled to the environment, is represented 
by the reduced density matrix (jTTJ) . In order to examine the reduced dynamics of the 
QIP, it is instructive to use the Kraus-map representation [51] . If the composite system 
was initially (i.e., at time i = 0) in the direct-product state, 

2" 

p(0) = p q (0) ® p cnv (0) = p q (0) ® 0»H ( 39 ) 
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then the reduced dynamics of the QIP has the following form (known as the Kraus map 

my- 

2 n 

Pq (i) = <% q (0)] = J2 K w ,{t)p^)KlAt), (40) 

u,v'=l 

where the Kraus operators K vv i (i) G M 2 m (C) are given by 

2 m 

K vv ,{t) = J2 U «' MlWl ( 41a ) 
i,i'=i vv ' 

2 n 

E Kl,{t)K„ v ,{t) = h m . (41b) 
i/,i/'=i 

It is well known (51] that there exist infinitely many different sets of Kraus operators, 
{Ki, . . . , K p } (where p G N is the number of operators in the set), that represent the 
same map $ (i.e., they evolve p q (0) in exactly the same way). Moreover, any Kraus map 
for an /c-level quantum system can be represented by a set of p < k 2 Kraus operators. 
That is, if the map is represented by a set of p' > k 2 Kraus operators, there always exists 
another representation with not more than k 2 operators. Therefore, for our system of m 
qubits and n environmental particles, the set of 2 2n Kraus operators ( 14 lap can always be 
transformed into another set of not more than 2 2m operators, representing the same map 
$. However, since we numerically study Kraus operators only for the case of n = m — 1, 
there is no practical need for such a transformation. 

In calculations, we use p(0) = |\E f )( x ^o| with |^ ) °f COB - F° r one qubit coupled 
to a one-particle environment, we use the notation \v = 1) = |+) and \v = 2) = |— ) 
and find K 12 (t) = K 22 (t) = and Kl^Kuit) + Kl x {t)K 2X {t) = J 2 . It is therefore 
sufficient to explore either Kn(t) or K 2 i(t). By evaluating the Kraus operators we can 
quantify the non-unitarity of the qubit dynamics. It is important to note that the non- 
unitary evolution is not only responsible for decoherence, but is also required to steer 
the information flow back to the QIP. The control field that restores coherence to the 
QIP necessarily employs the interaction with the environment and the corresponding 
non-unitary dynamics. We examine the time behavior of the Frobenius norm of 
the Kraus operator, ||i^ 21 (t) || Fr , that serves as a measure of non-unitarity. Figure [6] 
shows ||i^2i(^)||pr for both controlled and uncontrolled dynamics. In comparison to the 
uncontrolled evolution, the optimal control dramatically decreases the non-unitarity 
of the qubit dynamics during the gate operation, culminating in almost complete 
decoupling at the final time tf. We also see that, under the optimal control, ||i^2i(*)||pr is 
approximately symmetric about t ~ if/2. Inspecting the time derivative of the entropy, 
dS'vN/dt, we find that ||i^2i(i)||pr reaches the maximum at approximately the same time 
(just prior to if/2) when the fastest decrease in the qubit's entropy is observed, indicating 
the maximum flow of information into the QIP. 
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Table 1. The performance of the optimally controlled one-qubit Hadamard gate in 
the presence of various n-particle environments (7 = 0.02). Here, A max , tf, £, F, and 
S v n{U) are the maximum field amplitude, control duration, field fluence, gate fidelity, 
and final-time entropy, respectively. F 1= q denotes the gate fidelity obtained when the 
control field optimized for 7 = is applied to the system with 7 = 0.02. The initial 
state for the entropy computation is |Wq) of (fT9"l) . 



n 


1 


2 


4 


6 


A 


2.0 


4.0 


4.0 


2.5 


u 


25.0 


15.4 


25.0 


25.0 


£ 


20.0 


49.0 


55.5 


34.0 


F 


0.9995 


0.9975 


0.9935 


0.9786 


^7=0 


0.9063 


0.8829 


0.8133 


0.7723 


S v N(t{) 


9.0 x f0~ 8 


4.4 x 10~ 5 


4.7 x 10~ 4 


3.0 x 10~ 3 



6.3. One qubit coupled to a multiparticle environment 

We explore the performance of optimally controlled one-qubit gates in the presence of 
multiparticle environments described in section [2j Table [1] reports optimal control field 
parameters, fidelity, and final-time entropy for the one-qubit Hadamard gate coupled to 
n-particle environments (m = 1, n = 1,2,4,6, and 7 = 0.02). For n = 4, the values in 
Table [T] were obtained with the coupling scheme modeled by a two-dimensional lattice 
of flZJ); however, very similar results were obtained with the linear nearest-neighbor 
coupling scheme of (JSJ). 

The results obtained for n > 2 further illustrate the benefits of optimal controls 
which explicitly take into account coupling to the environment. The entropy dynamics 
indicate that for multiparticle environments the control employs the same mechanism 
of an induced coherence revival, as described above for n — 1. Fast and intense control 
fields significantly suppress the qubit-environment interaction during the gate operation 
and try to recover as much of the lost information as possible before the end of the 
control pulse. However, as the complexity of the composite system increases, it becomes 
more difficult to induce an almost perfect revival; therefore, the gate fidelity and final- 
time coherence decrease as n increases. This observation supports the conclusion that 
shorter-time controls (which do not rely on revivals) will be useful for environments with 
dense spectra. 

6.4. Two qubits with a one-particle environment 

For the QIP consisting of two qubits (m = 2), the target gate is CNOT of (I24bl) . The 
coupling constant between the two qubits is 712 = 0.1, while the coupling constant 
between each qubit and the single environmental particle (n = 1) is 713 = 723 = 7. 
Frequencies of the two qubits are uj\ = 1 and U2 = vr — 2.05 ~ 1.09519, and the 
frequency of the environmental particle is u 3 = (71 — 2.14)" 1 w 0.99841. The optimal 
control fields obtained for 7 = and 7 = 0.01 (shown in figure [7]) generate the CNOT 
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gate with fidelities of 0.9999 and 0.9798, respectively. When 7 = 0.01, the entropy for 
the uncontrolled evolution increases monotonically until t « 125 (reaching a maximum 
of approximately 0.6), whereas the optimal control field results in a much lower entropy, 
shown in sub-plot (b) of figure [HJ The same pattern of a partial revival at an intermediate 
time followed by an almost complete revival at t — if, seen for the one-qubit gates in 
figure is also present for the two-qubit gate, but on a longer time scale. For the 
CNOT gate's final-time coherence revival we find S v ^(t{) ~ 1.5 x 10~ 3 at tf = 121.1. 

We observe that the fidelity of the optimally controlled quantum gates decreases 
with increases in n (the number of environmental particles) and, even more significantly, 
m (the number of qubits in the QIP). This behavior arises due to the difference between 
the perfect control solution and an actual control field found by the optimization 
algorithm. According to an analysis of the control landscape for unitary transformations 
[52| 153] , the pernicious effect of control inaccuracies on the gate fidelity rapidly increases 
with the size of the system. If instead of the perfect control solution C*o(t), the actual 
field is Coif) + 5C(t), then instead of the perfect fidelity F — 1, one will obtain 
F — 1 — SF, where 5F oc 2 m \\5C(t) || 2 (here, || • || denotes an appropriate functional 
norm). As the number of interacting qubits, m, increases, the factor 2 m becomes more 
important. Moreover, as the complexity of the composite system increases (more qubits 
and/or environmental particles), the control error ||5C(£)|| will increase as well, as it 
will become more difficult to find a field that is very close to the perfect one. 

6.5. Can the state purity measure the gate fidelity? 

We found that obtaining a very high gate fidelity requires an almost complete coherence 
revival characterized by a very low final-time entropy. Is it then possible to rely on a 
characteristic of coherence (e.g., the final-time entropy or purity of the QIP state) as a 
measure of the gate quality, instead of measuring the distance between the actual and 
target gate transformations? The answer is definitely "no" because the restoration of 
coherence is a necessary, but not sufficient, condition for a high gate fidelity. There exist 
an infinite number of unitary, or almost unitary, transformations which nevertheless are 
very far from the target one. 

In order to further emphasize this point, we generalize the notion of the gate fidelity 
(as measured by the distance between the actual evolution operator U(t) and target 
transformation G) to all times < t < tf. Figure M shows this time-dependent fidelity 
F(t) and the entropy S v ^(t) for the optimally controlled two-qubit CNOT gate (with 
7 = 0.01 and tf = 121.1). We see that the minimum of the entropy occurs at a time 
£s min « 119 (i.e., before tf) when F(t) is still quite low, and that at the time interval 
between tg mln and tf, while the fidelity F(t) rapidly increases to achieve its final-time 
value F ~ 0.9798, the entropy slightly increases as well. This example shows that 
fidelity and coherence do not always correlate and that a very low value of the entropy 
does not always result in a correspondingly high value of the gate fidelity. According 
to this analysis, a strategy of maximizing the state purity [381 SO] does not ensure the 
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generation of target quantum gates with the highest possible fidelity. 

7. Robustness of optimally controlled gates to system variations 

We observed that applying the control field optimized for the closed system (7 = 0) 
to the coupled one (7 = 0.02) results in a significant decrease in the gate fidelity. 
Analogously, we find that applying the control field optimized for the case of a one- 
particle environment (n = 1) to systems with n > 2 environmental particles also has 
a strong detrimental effect on the gate fidelity. These results are part of a broader 
analysis of the robustness of optimally controlled quantum gates to different types of 
system variations. 

We address some aspects of this issue by considering the one-qubit Hadamard gate, 
with a fixed number n of environmental particles (n = 1,2,4), and finding an optimal 
control field for a specified set of system parameters: the coupling constants 7^ given 
by © (with 7 = 0.02) and frequencies u>i given by (122]) for n = 1 and (1231) for n > 2. 
Then we apply this control field to an ensemble of systems with normal variations in 
either coupling constants 7^ or frequencies o>j and analyze how the uncertainties in the 
system parameters affect the gate fidelity F and final-time entropy SVn^)- Although 
the dependence of F and S v ^(t{) on the coupling constants and frequencies is highly 
non-linear (which implies that the distributions of F and S^N^f) will not be normal), 
our statistical analysis employs only mean values and standard deviations, given by 
F = L- 1 J2r=i F r and a F = [L~ l Y%=i( F r - F) 2 } 1 ' 2 , respectively, for the gate fidelity 
F, and similarly for the final-time entropy S v ^(t{). The summation is over all elements 
of the ensemble (ensemble sizes L of the order of 10 5 are used in the calculations). 

7.1. Variation of the coupling constants 

The value of each non-zero coupling constant 7^ (given by (jSJ) with 7 = 0.02) is 
individually replaced by a value randomly selected from a normal distribution with 
a mean 7 = 0.02 and a standard deviation 07 = 7/8 = 0.0025. The statistical analysis 
of the corresponding distributions of the fidelity and final-time entropy is reported in 
table El and frequency histograms of these distributions are shown in figure [9j These 
results demonstrate a high degree of robustness of the performance of the optimally 
controlled gate to relatively large variations in the strength of the system-environment 
coupling. On average, there is practically no decrease in the fidelity and entropy, and 
the relative width of the fidelity distribution, crp/F, is by several orders of magnitude 
smaller than 0^/7. Interestingly, if the control field optimized for 7 = 0.02 is applied to 
the closed system with 7 = 0, this results in a relatively high fidelity (e.g., F = 0.9989 
for n=l). The standard deviation o> rises with the increase in the number of 
environmental particles. We also see that the distributions of F and SVn^O are more 
symmetric for n = 4 than for n = 1. 
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Table 2. Fidelity and entropy data for the one-qubit Hadamard gate applied to 
an ensemble of systems with normal variations in the coupling constants 7^ and 
frequencies Wj. Columns of F and S v yi(t{) contain fidelity and final-time entropy 
values, respectively, for the original system parameters: 7 = 0.02 and frequencies 
given by ([2"2")) for n — 1 and (|2"3")l for n > 2. Columns of F and SVn contain mean 
values of fidelity and final-time entropy, respectively, over the ensemble, while <jf and 
<7s vN are the respective standard deviations. 
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0.9896 
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0.9935 


0.9884 


4.5 x 10~ 3 


4.7 x 10- 


-4 


1.7 x 10 


-3 


1.8 x 10" 


-3 



7.2. Variation of the frequencies 

The value of each frequency u>i (given by ff22l) forn = 1 and ff23l) for n > 2) is individually 
replaced by a value randomly selected from a normal distribution with a mean uJi = uoi 
and a standard deviation a Ui = Ui/25. The statistical analysis of the corresponding 
distributions of the fidelity and final-time entropy is reported in table [2j and frequency 
histograms of these distributions are shown in figure [TDl It is well known [2"9"l [5UI [52] 
that a high degree of quantum control may be achieved through the complex interference 
of evolution pathways. This interference strongly depends on the relative phases of all 
pathways, and these phases in turn depend on the transition frequencies of the system. 
Therefore, we would expect the optimal gate performance to be much more sensitive 
to variations in the frequencies than to changes in the coupling constants. The results 
presented in table [2] and figure [10] corroborate this expectation. Still, the robustness 
of the optimal gate performance to frequency fluctuations is tolerable. Moreover, the 
degree of robustness for systems with two and more environmental particles (n > 2) is 
even higher than for n — 1. 

8. Conclusions 

This work demonstrates the importance of OCT in designing quantum gates for use in 
QC, especially in the presence of a decohering environment. The model studied here 
represents a realistic system of interacting qubits and is relevant for various physical 
implementations of QC. High quality optimal solutions obtained in the presence of 
unwanted couplings also exhibit a significant degree of robustness to random variations 
in the system parameters. The analysis of the system dynamics reveals control 
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mechanisms which employ fast and intense time-dependent fields to effectively suppress 
the qubit-environment interaction via dynamic shifting of the energy levels and achieve 
an almost full coherence recovery via an induced revival. 

The results reported in this paper further support the use in QC applications 
of laboratory closed-loop optimal controls employing learning algorithms and intense 
ultrafast fields [29], [301]. In the area of molecular dynamics, the utility of optimal 
control methods was first demonstrated theoretically in very simple model systems; 
nevertheless, these methods were later applied with great success in the laboratory to 
complex molecules [29J. Similarly, we expect that the optimal control of quantum gates, 
the usefulness of which was demonstrated here for a relatively simple environment model, 
will be also effective for real quantum information systems. A successful application 
of optimal control methods to the generation of high-fidelity quantum gates in the 
laboratory will be an important step towards achieving error thresholds required for 
fault tolerant QC |S1E]- 

This work may be further advanced with the use of the control-mechanism analysis 
[54] to explore the detailed dynamics of the decoherence management process in 
optimally controlled quantum gates. Methods of landscape analysis [521 [531 EE] may 
be employed to investigate how optimal controls are deduced and study the effects of 
control errors in the context of non-unitary dynamics of open quantum systems. 
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Figure 1. The time-evolution of the entropy S v ^(t) for the uncontrolled system of one 
qubit coupled to a one-particle environment, with 7 = 0.02, toi = 1, and various values 
ofw2- Solid lines: LO2 — (ir — x)^ 1 ; dashed lines: LO2 = K — x (with x = 2, 2.1, 2.14). The 
initial state is I'J'o) of (fT9|) . For a given value of 7, closer frequencies cui and LO2 enhance 
the interaction between the qubit and environment, causing stronger decoherence and 
longer revival times. 




Figure 2. The time-evolution of the entropy S v ^(t) for the uncontrolled systems of 
one qubit coupled to n-particle environments: n = 2 (solid line), n = 4 (dashed line), 
and n = 6 (dotted line). The coupling constant is 7 = 0.02. Frequencies of the qubit, 
u>x, and the environmental particles, ujj (j — 2, . . . ,n + 1), are given by (|23p . The 
initial state is |* ) of (fH)]) . 
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Figure 3. The gate fidelity F versus the coupling constant 7, for optimally controlled 
one-qubit gates: Hadamard (solid line), identity (dashed line), and phase (dotted line). 
Each one-qubit gate is coupled to a one-particle environment. Values of 7 range from 
to 0.02 in increments of 0.001. 
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Figure 4. Optimal control fields C(t) versus time, for one-qubit gates: (a) Hadamard, 
(b) identity, and (c) phase. Each one-qubit gate is coupled to a one-particle 
environment (7 = 0.02). 
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Figure 5. The von Neumann entropy S v -^(t) versus time, for optimally controlled 
one-qubit gates: (a) Hadamard, (b) identity, and (c) phase. Each one-qubit gate is 
coupled to a one-particle environment (7 = 0.02). The initial state is \^>q) of (p~9|) . 
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Figure 6. The time-evolution of the Kraus operator's norm, ||-K2i(t)||Fr> for the system 
of one qubit and one environmental particle: uncontrolled evolution with 7 = 0.1 
(dotted line), uncontrolled evolution with 7 = 0.02 (dashed lines), and controlled 
evolution, under the optimal control field generating the Hadamard gate, with 7 = 0.02 
(solid line). The initial state is |* ) of CUD - 
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Figure 7. Optimal control fields C(t) versus time, for the two-qubit CNOT gate with 
7 = 0.01 (solid line) and 7 = (dotted line). 
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Time 

Figure 8. The time-evolution of (a) fidelity F(t) and (b) von Neumann entropy S v ^(t) 
during the optimally controlled CNOT gate operation (m = 2, n = 1, and 7 = 0.01). 
These results demonstrate that a high degree of coherence (quantified by the entropy) 
does not ensure a correspondingly high value of the gate fidelity. The initial state for 
the entropy computation is |^o) of (fT9|) . 
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Figure 9. Frequency histograms for the gate fidelity and final-time entropy 
distributions, obtained when the control field optimized for the Hadamard gate with 
7 = 0.02 is applied to an ensemble of systems with normal variations in the coupling 
constants 7y. The distribution for each non-zero 7^ is normal with a mean 7 = 0.02 
and a standard deviation a 1 = 7/8 = 0.0025. Sub-plots include frequency histograms 
of (a) the fidelity distribution for n — 1, (b) the entropy distribution for n = 1, (c) 
the fidelity distribution for n = 4, and (d) the entropy distribution for n = 4. Note 
the axes scale differences in the sub-plots. Table [2] reports statistical data for these 
distributions. 
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Figure 10. Frequency histograms for the gate fidelity and final-time entropy 
distributions, obtained when the control field optimized for the Hadamard gate with 
transition frequencies o>j given by (|22[) for n = 1 and ([23]) for n > 2 is applied to 
an ensemble of systems with normal variations in the transition frequencies. The 
distribution for each transition frequency is normal with a mean = u>i and a standard 
deviation a UJi — Wj/25. Sub-plots include histograms of (a) the fidelity distribution for 
n = 1, (b) the entropy distribution for n = 1, (c) the fidelity distribution for n — A, and 
(d) the entropy distribution for n = 4. Note the axes scale differences in the sub-plots. 
Table [2] reports statistical data for these distributions. 



